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Abstract. We investigate the radio spectral energy distributions (SEDs) of young star-forming galaxies and how 
they evolve with time. The duration and luminosity of the nonthermal radio emission from supernova remnants 
(SNRs) are constrained by using the observational radio SEDs of SBS 0335—052 and I Zw 18, which are the two 
lowest-metallicity blue compact dwarf galaxies in the nearby universe. The typical radio "fluence" for SNRs in 
SBS 0335—052, that is the radio energy emitted per SNR over its radiative lifetime, is estimated to be ~ 6- 
22 X 10^^ W Hz"^ yr at 5 GHz. On the other hand, the radio fluence in IZwl8 is ~ 1-3 x 10^^ W Hz"^ yr at 
5 GHz. We discuss the origin of this variation and propose scaling relations between synchrotron luminosity and 
gas density. We have also predicted the time dependence of the radio spectral index and of the spectrum itself, for 
both the "active" (SBS 0335—052) and "passive" (I Zw 18) cases. These models enable us to roughly age date and 
classify radio spectra of star-forming galaxies into active/passive classes. Implications for high-z galaxy evolution 
are also discussed. 
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1. Introduction 

Radio emission from galaxies is known to be connected 
with star formation activity. Two radiative processes are 
responsible: thermal free-free radiation from ionized gas 
in H II regions and synchrotron radiation from relativistic 
electrons spiraling in magnetic fields (Condon 11992(1 . The 
former originates from ionized gas around massive stars, 
while the latter comes from supernova remnants (SNRs) 
whose progenitors are massive stars. Since massive stars 
have short lifetimes, both thermal and nonthermal emis- 
sion should trace the current star formation rate (SFR). 

Nevertheless, the physical basis of the radio emission 
is not fully established. The observed radio flux in lumi- 
nous evolved starbursts and normal spiral galaxies is dom- 
inated by nonthermal emission, but not from the short 
timescale radiation of discrete SNRs; more than 90% of 
it is probably due to the diffusion of cosmic ray electrons 
in the galaxy disk (Condon 11992^ over timescales of 10^- 
10* yr (Helou & Bicav [TM5|l . Indeed, if the theoretically 
estimated adiabatic timescale (e.g., Woltier ri972|l and the 
observationally obtained T,-D (surface brightness vs. di- 
ameter) relation for discrete SNRs (e.g., Clark & Caswell 
I1976|l are used to derive the nonthermal emission in galax- 



ies (e.g., Biermann 119761 Ulvest ad 11982(1 . the emission is 
underpredicted by a factor of 10 or so. 

In very young galaxies with ages ^ 10 Myr, however, 
the cosmic-ray diffusion mechanisms have not yet had 
a chance to dominate the radio emission. Nonthermal 
synchrotron radiation from discrete SNRs would be ex- 
pected to dominate over the diffuse component. Moreover, 
in galaxies with starbursts younger than ^ 3 Myr, there 
should be very little synchrotron emission since the SNRs 
have not yet exploded. Indeed, some galaxies with a deficit 
of nonthermal emission are observed (Roussel et al. 12003(1 
and appear to be due to starbursts observed within a few 
Myrs of their onset, after a long period 100 Myr) of 
quiescence. 

In this paper, we investigate theoretically the time evo- 
lution of the radio spectral energy distribution (SED) in 
young star-forming galaxies. In principle, this is a straight- 
forward process because we know that the radio SED can 
be determined by the star formation history, for which 
we have already developed the formalism (Hirashita et al. 
IT002.: Hirashita & Hunt IMHjl . In practice, it is a diffi- 
cult exercise because of the unknown nature of the time 
evolution of the nonthermal radio component. Hence, we 
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adopt nearby blue compact dwarf galaxies (BCDs) as an 
observational sample to constrain our theoretical model 
of radio emission. Most BCDs have a young age of cur- 
rent star formation (although most of them had previous 
starburst episodes), and a low mctallicity (e.g., Izotov & 
Thuan p004 ) . This implies that they are relatively un- 
evolved chemically and may provide a reasonable tem- 
plate to investigate star formation properties of young 
galaxies. The BCDs as a class have radio spectra with 
different properties than those of evolved luminous star- 
forming galaxies (Klein et al. ,.1991j : flatter radio spectra of 
BCDs indicate that they have a lower fraction of nonther- 
mal emission than more luminous systems. Hence, they 
may lack the diffuse emission that characterizes the radio 
spectrum of larger disk galaxies, and thus are ideal targets 
for constraining our models. 

Two classes of star-formation activity in BCDs have 
recently emerged observationally, as proposed by Hunt 
et al. H2003a|l and Hunt & Hhashita (2006, m prepara- 
tion). They argue that the star-formation modes in the 
two most metal-poor galaxies, SBS 0335—052 and I Zw 
18, are very different, in spite of their similar metallicities 
(12 + log(0/H) = 7.3 and 12 + log(0/H) = 7.2, re spec- 
tively; Skillman & Kennicutt 1993; Izotov et al. lTM^ . The 
major star-forming region of SBS 0335— 052 is compact 
and dense (radius rsF ^ 40 pc, number density n <^ 600 
cm-3; Dale et al. lMIIl Izotov & ThuanUHMIl- Moreover, 
SBS 0335— 052 hosts several super star clusters (SSCs), 
detectable H2 emission lines in the near-infrared (NIR) 
(Vanzi et al. l2()()()(l . a large dust extinction {Ay ^ 16 mag; 
Thuan et al. [T^M Hunt et al. IWiTl Plante & Sauvage 
I20()2|l . and high dust temperature (Hunt et al. 2001, Dale 
et al. IMJTl Takeuchi et al. IMT5|l . On the contrary, the 
star-forming regions in IZwl8 are diffuse (rsp ^ 100 pc, 
n ^ 100 cm^'^), and contain no SSCs. NIR H2 emission 
has not been detected (Hunt et al., private communica- 
tion), and the dust extinction is moderate (Ay ~ 0.2 mag; 
Cannon et al. 2002). We consider a region with such prop- 
erties as "passive" following Hunt et al. H2003a|l . while 
a SBS 0335— 052-like star-forming property is considered 
"active". The similar metallicities of SBS 0335— 052 (ac- 
tive) and I Zw 18 (passive) imply that the chemical abun- 
dance is not a primary factor in determining the star- 
forming properties. We argue that the compactness of 
star-forming regions, which affect gas density, gas dynam- 
ics, and so on, is important in the dichotomy of active and 
passive modes. 

In this paper, we extend this active-passive classifica- 
tion into the radio regime by focusing on SBS 0335— 052 
and IZwl8 as representatives of metal-poor BCDs, and 
also of two star-forming modes in BCDs: "active" and 
"passive". This paper is organized as follows. First, in 
Sect. 12 we explain the model that describes the evolu- 
tion of radio SEDs. Some basic results are given in Sect.O 
where observational data of SBS 0335—052 and I Zw 18 are 
compared with the results for active and passive modes. 
The constraints obtained here are compared with theoret- 
ical descriptions of the physical processes of synchrotron 



radiation in Sect.01 In Sect. [S] we generalize the results of 
our models with predictions for the time evolution of radio 
SEDs distinguishing between active and passive modes, 
and discuss implications for high-redshift star formation. 
Finally we give our conclusions in Sect. El 



2. Model description 

The radio emission from galaxies is interpreted to be com- 
posed of two components fCondon ll992|l : thermal free-free 
(bremsstrahlung) emission and nonthermal synchrotron 
radiation. The former requires treatment of ionized re- 
gions around star-forming regions and the latter is related 
to SNRs. Our models include these processes as described 
below. We approximate the dominant star-forming region 
in a BCD to a single zone with spherical symmetry. 



2.1. Thermal component 



The free-free radiation in ionized regions is responsible for 
the flat thermal component of radio SEDs. The luminosity 
of the thermal component is proportional to the number of 
ionizing photons (with energy larger than 13.6 eV) emitted 
per unit time, A^ion- In order to calculate iVion, we need the 
stellar ionizing photon luminosity and the stellar lifetime. 

We use the fitting formulae of Schaerer ^2002) for the 
stellar lifetime (denoted as t^, in Schaerer i2002 ) and the 
number of ionizing photons emitted per unit time Q{m) 
(see Table 6 of Schaerer I2002|l . Those quantities are cal- 
culated as functions of stellar mass at the zero age main 
sequence, m. We adopt Z = (zero metallicity) for the 
stellar properties because we treat a metal-poor phase. If 
we adopt the solar metallicity instead, A'ion is 2 times 
smaller (Schaerer 2002). Considering those two extreme 
metallicities, we consider that the uncertainty caused by 
the stellar metallicity is within a factor of 2. 

The evolution of A^ion as a function of time t is calcu- 
lated by 



Nionit) 



dm I dt' Q{m)(j){m)ip{t ~t') 
'0 



(1) 



where ip{t) is the SFR at t, and (j){m) is the initial mass 
function (IMF; the definition of the IMF is the same as 
that in Tinslev ri98()|> . In this paper, we assume a Salpeter 
(Salpeter 11955(1 with stellar mass 



TO oc m 



-2.35 



IMF 

range of 0.1-100 M©. With this IMF, the mass fraction 
of stars more massive than 5 Mq is 0.18. The calculated 
radio luminosities are roughly scaled with this fraction of 
massive stars if another IMF is adopted. The functional 
form of the SFR ?/;(t) is specified in Sect. 12.31 

The number of ionizing photons can be related to the 
thermal radio luminosity at the frequency Ljjj(i') as 
(Hunt et al.EOnil hereafter H04; valid for 10000 K T 
20000 K) 



1.58 X 10 
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T 



W K 



-0.31 



1052 s- 



W Hz" 



(2) 



where n(H"'") and n(He^) are the number densities of ion- 
ized hydrogen and ionized heHum, respectively (we assume 
n(He+) = 0.08n(H+) in this paper; H04), and T is the gas 
temperature. For the gas temperature in the ionized re- 
gion, we assume a typical value of metal-poor H ii regions: 
T ~ 2 X lO'^ K (SkiUman & Kennicutt [T99J Izotov et al. 
11999 ). The change of the temperature from 2 x 10* K to 
1 X 10** K has only a minor effect on L\ and causes a fac- 
tor of 2.5 difference in the free-free optical depth estimated 
below (Eq.OJ. 

A part of the radio emission is self-absorbed by the 
ionized gas. We estimate the optical depth of the free-free 
radiation, as 



Tff ~ 0.328 



104 K 



-1.35 



EM 



10^ pc cm~ 



GHz 



-2.1 



(3) 



(H04), where the emission measure EM is estimated by 
assuming a constant electron number density n^: EM = 
n\l {£ is the length of the ionized region along the line of 
sight). The length £ is estimated by dividing the volume 
with the surface of the ionized region. We thus obtain the 
emission measure as 



EM = -nir^ 



(4) 



where n is the radius of ionized region defined in Sect. 12.41 
In the ionized region, the following condition is assumed 
to be satisfied for the charge neutrality (H04): 



He = n(H 



?i(Hc^ 



l.OSnn 



(5) 



where rifj is the number of hydrogen nuclei. Since the ion- 
ized gas is responsible for both emission and absorption, 
it is reasonable to assume that the absorbing medium 
is intermixed with the emitting material. In this case, 
the monochromatic luminosity of the thermal component 
emerging from the ionized region becomes 



r / ^ 1 -CXp(-Tff) 



Tff 



2.2. Nonthermal component 



Nonthermal radio emission from galaxies originates from 
SNRs. Since we consider young star-forming regions, we 
only treat Type II supernovae (SNe II), whose progenitors 
have short lifetimes^. The rate of SNe II as a function of 
time, 7(i), is given by 



7(i) = / ^{t 

J 8 Mq 



(m) dm . 



(7) 



First we use a "standard" estimate for the intrinsic 
nonthermal radio luminosity at frequency L[Jjg(z^) by 

^ Although see Mannucci et al. II2005II for the idea that Type 
la supernovae are important in young systems. 



relating it to 7(t) (e.g., Condon [T992|l 

-1/17 

LlU:.) = 1.2 X 10^ 



j^O^i erg 



1 cm 



^ TwHz^^ 

1 yr-iy V408 MHz^ 



(8) 
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where E'en is the explosion energy of a SN II (E'sn = 10 
erg is assumed in this paper), and a is the spectral index. 
We designate this model for the nonthermal component 
as the "s-model" ("s" stands for standard). However, this 
formula is obtained by assuming the S — I? (surface bright- 
ness vs. diameter) relation derived from the Galactic SNRs 
(Clark & CasweU 1976; Ulvestad 1982) and the adiabatic 
lifetime Tad as time during which the remnant emits in 
the radio. The adiabatic timescale is estimated as (Woltjer 
^372, Condon J 992): 



Tad =^ 3.4 X 10* 



105 



\ 4/17 

1 erg J \ 1 



riH 



cm 



-9/17 



y. (9) 



The SN II rate inferred from the nonthermal radio lumi- 
nosities by using Eq. © is generally too high (Condon 
119921 and references therein; see also Sect. In other 
words, the radio luminosity given by the dependence on 7 
in Eq. ||HJ) is an underestimate. Thus, we need to reexam- 
ine the underlying physical assumptions. Condon l)1992|l 
suggests that the discrepancy arises because SNRs pro- 
duce more than 90% of the nonthermal radio emission 
long after the individual SNRs have faded out and the 
electrons have diffused throughout the galaxy. However, 
this may not be true in young starbursts, hence warrant- 
ing an investigation of the possibility that individual SNRs 
have a radio lifetime longer than the adiabatic timescale 
(Ilovaisky & Leaueux ll972|l . 

In order to constrain the radio energy emitted by 
a SNR over its entire lifetime - the "fluence" or time- 
luminosity integral - we treat the duration (rnt) and the 
luminosity at a frequency z/q (^nt) as parameters to be de- 
termined from observational constraints. Then, the non- 
thermal luminosity is written as 



(6) Ll,(y) = Ur^ai^/v^r ■ 



(10) 



This model of the nonthermal component is called "p- 
niodel" ("p" stands for parameterization). We take vq = 
5 GHz unless otherwise stated. We adopt a typical spectral 
index of the Galactic SNRs: a ~ —0.5 (Clark & Caswell 
ITOTSj) . Hunt et al. H2005|l mention that a = -0.4 and 
—0.5 fit the radio spectrum of IZwl8 but that the fit is 
significantly degraded for a — —0.8. 

After accounting for the absorption, the nonthermal 
radio component emitted by the star-forming region is es- 
timated as 



int(i')=£;(rff)i" 



(11) 



where E{Tg) is the escaping fraction of the radio photons, 
and i is the label indicating "s" for the s-model and "p" for 
the p-model. We consider two cases for E{Tg) depending 
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on the geometry of the absorbing material relative to the 
radiation source. One is the "screen geometry" , where the 
absorbing material surrounds the radiation source, and 
the other is the "mixed geometry" , where the source and 
the absorber are mixed homogeneously: 



Escirtf) = exp(-Tff) 
E, 



(screen) , 



1 - cxp(-rff ) . 
(Tff) = (mixed). 



(12) 



2.3. Star formation rate 

Stars form as a result of the gravitational collapse of a 
gas cloud. Therefore, it is physically reasonable to relate 
the SFR with the free-fall (or dynamical) timescale of gas 
(Elmegreen BUPDI - We consider a star-forming region with 
an initial number density of neutral hydrogen tiho- The 
free-fall time, tg, is estimated as 

/ Stt 4.35 X 10^ 
iff ^a/ttttt^ yr, (13) 



32Gp 



where G is the gravitational constant, and p is the mass 
density of the gas. In Eq. (|13|l . we have used p = /Ltmn'^HO) 
where ttih is the mass of a hydrogen atom, and the factor 
H is for the correction for helium. We assume that p, — 1.4 
in this paper. 

The SFR, ip{t), basically scales with the gas mass di- 
vided by the free-fall time: 

espMgas 



--fit), 



(14) 



where esF is the star formation efficiency defined as the 
conversion efficiency of a gas into stars over a free-fall time, 
Afgas is the total gas mass available in the star-forming re- 
gion, and the nondimensional function f{t) specifies the 
functional form of the SFR. Numerically, the SFR is esti- 
mated as 



V'(i) ^ 0.230 



\o.iJ \ w Mq 

X fit) Mq yr-1 . 



;)( 



100 cm- 



1/2 



(15) 



We set the zero point of time t at the onset of star for- 
mation in the star-forming region. We define the termina- 
tion time of star formation as ti such that ri{ti) = rsF, 
where tsf is the radius of the star- forming region (i.e., tf 
is the time when the entire star-forming region is ionized; 
see Sect. 12.41 for the definition of n). If the entire region 
is never ionized, ti is taken to be infinity. We note that 
tff < ti is satisfied for the conditions adopted in this pa- 
per. Then, we investigate the following functional form for 

fit)-- 

fit) = fcit) ^ 



exp(— esF^/iff) if < t < ^i, 
if<<Oort>ti, 



(16) 



which is called "continuous SFR" . In order to represent an 
extreme starburst, we also adopt another functional form 
for the SFR: 



fit) = /b(i) 



s if < t < tff, 

if < < or i > tff, 



(17) 



where the constant s is adjusted to satisfy 



1 



= — / exp(-esFiAfT) 
tff Jq 

^ [1 - exp(-esFii/tfr)] , 

esF 



(18) 



where ti is determined by the continuous SFR model. 
This assures that the total stellar mass formed in a star- 
forming region is the same between the two models. In the 
burst model, all the star formation occurs on a dynamical 
timescale. For the star formation efficiency, we assign a 
value esF = 0.1 (e.g., Inoue et al. l2UUU|l . 

The initial hydrogen number density can be related to 
the gas mass in the star-forming region as 



47r 



(19) 



Then, the numerical value of the number density is esti- 
mated as 



an f '^SF \ ( A/gas \ _3 

nno ^ 69 ( —. | ( . I cm 



VlOOpcy Mq 



2.4. Radius of ionized region 



(20) 



In order to treat the range of evolution variations from 
deeply embedded Hll regions to normal Hli regions, it 
is crucial to include pressure-driven expansion of H ii re- 
gions. To focus on the radio SED, we adopt a simple an- 
alytical approximation, which divides the evolution of an 
Hii region into two stages: the first stage is the growth of 
ionizing front due to the increase of ionizing photons, and 
the second is the pressure-driven expansion of ionized gas. 
The expansion speed of ionizing front in the first stage is 
simply estimated by the increasing rate of the Stromgren 
radius. The Stromgren radius rso under the initial density 
is defined by the following relation: 

' Nion, (21) 



-rlonoriHoQ!*^^ 



where a(2) is the recombination coefficient excluding cap- 
tures to the ground (n = 1) level, and tt-o is estimated by 
Eq. jSJ with riH — ^^Ho• The increase of rgo is caused by 
the accumulation of ionizing stars. Roughly speaking, as 
long as rgo (the increase rate of the Stromgren radius) is 
larger than the sound speed of ionized gas, Cn (we as- 
sume Cii = 10 km in this paper), the ionizing front 
propagates before the system responds hydrodynamically. 
Therefore, we neglect the hydrodynamical expansion if 
rso > Cii, and assume that the radius of ionized region is 
equal to the Stromgren radius: n = rso. 

Once rs < Cu is satisfied, pressure-driven expansion is 
treated. Since the density evolves, we define the Stromgren 
radius rs under the current density: 



(22) 



47r n 

— rgrie'^Ha^ ' = A'ion ■ 

In this situation, the growth of the ionizing region is gov- 
erned by the pressure of ionized gas and the luminosity 
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change of the central stars has minor effects. Therefore, 
the following equation derived by assuming a constant lu- 
minosity fSpitzer I1978|l approximately holds: 



Table 1. Models of Radio SEDs 



(23) 



where pi and pu are the gas densities outside and inside 
of the ionized region, respectively. Here, pu/pi is approx- 
imately estimated by 



pu 
Pi 



3/2 



(24) 



where rjo is evaluated by the Stromgren radius under the 
current ionizing photon luminosity and the gas density of 
pi. Since nn is defined in the ionized region, we relate it 
with Pu as 



Pu = p.munu ■ 



(25) 



The density in the neutral region is assumed to be con- 
stant: 



Pi = ^J-munuo ■ 



(26) 



Equations (|23|) and (|24|l are numerically integrated to 
obtain ri as a function of t. For the convenience of in- 
terpretation, some approximate analytical solutions are 
described in Appendix IXI 

The dynamical expansion is treated as long as ri < rs . 
When the SFR declines significantly, rs begins to decrease. 
Thus, Ti > rs can be realized at a certain time. When 
Ti > rs, n is replaced with rs and nu is fixed; that is, we 
finish treating the dynamical expansion. 

3. Observational constraints 

Here we use observations to constrain some of the phys- 
ical parameters in our models. In particular, the radia- 
tive energy -time-luminosity integral- of the nonthermal 
synchrotron component is the most important parameter 
to be obtained in this paper. Two representative metal- 
poor BCDs, SBS 0335-052 and I Zw 18 (Hirashita & Hunt 
I2004|l . are used here. In our previous works (Hirashita et 
al.lSnD?; Hirashita & Hmit 200T), SBS 0335-052 has been 
used as a proxy for a genuinely young galaxy in which the 
mass of underlying old population is negligible (Vanzi et 
al. I2000|l . The same may be true for IZwlS in which the 
mass fraction due to an underlying evolved stellar popula- 
tion (^2 Gyr) is not more than - 20% (Hunt et al. ^003b). 
We adopt 54.3 Mpc for the distance of SBS 0335-052 
(Thuan et al. lTMTjl and 12.6 Mpc for I Zw 18 (Ostlin 2000). 

We calculate the time dependence of the radio SED, 
L{v), by summing the thermal and nonthermal compo- 
nents: 



L{l') = Lth(j^) + Lntiv) 



(27) 



Each of the following three items has two selections as 
summarized in Tabled 

(i) nonthermal emission of a SNR: standard / parameter- 
ized (Sect. 



Model 


rO a 


E{m) " 


SFR" 


A 


S 


s 


c 


B 


S 


s 


b 


C 


S 


m 


c 


D 


S 


m 


b 


E 


P 


s 


c 


F 


P 


s 


b 


G 


P 


m 


c 


H 


P 


m 


b 



See Section 12.21 The standard and parameterized models 
are denoted as "s" and "p" , respectively (Eqs. |H1 and 1101 
respectively). 

^ See Section The screen and mixed models are denoted 
as "s" and "m", respectively (Ea. ll2ll . 
See Section 12.31 The continuous and burst SFRs are de- 
noted as "c" and "b", respectively (Egs. lT^ and llTI respec- 
tively). 



(ii) optical depth of nonthermal emission: screen / mixed 
(Sect.E2J| 

(iii) SFR: continuous/burst fSect. 1231) 

In total we have 8 models, which are labeled as Models 
A-H as shown in Tabled 

3.1. SBS 0335- 052 

We adopt rsp = 20 pc and rino = 7 x 10^ cm~'^ (i.e., 
Mgas = 8.1 X 10^ Mq- ts = 5.2 X 10^ yr). We selected 
these values after surveying the full range of parameters 
so that the observable quantities (SFR, nu, and EM) are 
consistent with observations at the age estimated by Vanzi 
et al. H2000|l . The time evolution of the relevant quantities 
is shown below to verify whether those values indeed re- 
produce the observations. The radius of the star-forming 
region is also consistent with Takeuchi et al. (2005), who 
reproduce the observational FIR SED of SBS 0335-052. 

First, we evaluate the continuous SFR. The time evo- 
lution of n, nu, EM, and SFR is shown in Figures QJl, 
b, c, and d, respectively. The SFR averaged for 5Myr 
is ~ 0.8 Mq yr~^ with the total stellar mass formed of 
4 X 10^ Mq. The SFR of stars more massive than 5 Mq 
(SFR>5 Mq), which is less affected by the assumed IMF, 
is estimated to be 0.14 Mq yr~^, which is consistent with 
the estimate by H04 (SFR>5 Mq ^ 0.13-0.15 Mq yr'^). 
At t — 5 Myr, nu — 7 x 10^ cm~^, which is also com- 
patible with the observational results (Izotov et al. 119991 
H04). 

The radius of the ionized region ri monotonically in- 
creases until t = 3.7 Myr because of the pressure-driven 
expansion and the increase of iVion. At t — 3.7 Myr, the 
entire star-forming region is ionized. Until then, the in- 
creasing behavior of is approximated by rj oc t^^'^ as 
derived in Appendix IXI In Figure^, we show the slope 
of ri oc i^/^ (dotted line), which is roughly consistent with 
the calculated result (solid line). The initial evolution of 
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Fig. 1. Time evolution of basic quantities concerning radio emission, with the solid and dashed lines corresponding 
to the continuous and burst SFR, respectively. The initial conditions are selected to be the "active" mode, consistent 
with the observations of SBS 0335—052 at t ~ 5 Myr. a) Radius of ionized region n as a function of time, b) Number 
density of protons in the ionized region nn. c) Emission measure, d) Star formation rate. The dotted lines show the 
analytical approximation: n oc i^/^, tih oc t^^/*", and EM oc t~'^^'' . 



density is also explained by the analytical relation derived 
in Appendix as shown in Figure^ (nn oc As a 

result, the emission measure decreases as EM cx t^'^l'^ as 
shown in Figure^. 

If we adopt the burst SFR (s = 5.1; Eq.^J, we obtain 
the time evolutions as depicted in Figure^ As mentioned 
above, value of s is determined to ensure that the total 
stellar mass formed is the same between the burst and con- 
tinuous modes. Since the star formation stops on a short 
timescale, the driving force of the dynamical expansion 
decreases more rapidly than in the case of the continuous 
SFR. Thus, the maximum radius of the ionized region is 
smaller and the density is kept higher. This high density 
makes the EM higher in the burst SFR, and the EM de- 
rived observationally by H04 (EM = 3 - 7 x lO'^ pc cm~^) 
is more consistent with the burst SFR than with the con- 
tinuous SFR sXt ^ h Myr. Although the burst scenario is 
more consistent with the observations of SBS 0335— 052, 



in what follows we examine the two cases as extremes for 
the SFR. 

We show the radio SEDs calculated by various models 
listed in Tabled and examine whether or not the model 
predictions agree with observational data. The radio data 
are shown in Figure 12 the squares are from H04, and the 
cross is from Dale et al. (P7IlT|l . The SED at * = 5 Myr is 
presented in Figure 13 where the dotted and dashed lines 
represent the nonthermal and thermal components, re- 
spectively. The sum of those two components is shown by 
the solid line in each panel. The panel (a) shows the result 
of Model A, in which we adopt the "standard" model for 
the nonthermal radio emission. In this case, the thermal 
emission dominates the radio SED. The model prediction 
is inconsistent with the observational SED in the follow- 
ing two aspects: (i) the theoretical flux is systematically 
smaller than the observed one; (ii) the flat SED, domi- 
nated by the free- free component, is inconsistent with the 
observational SED at 4 GHz, which requires a promi- 
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Table 2. Best Fit Parameters for the Nonthermal 
Component (SBS 0335-052) 



Model 


;ntrnt(5 GHz) 


x' 


xi 




(W Hz-^ yr) 






E 


7.6 X 10^^ 


6.8 


1.4 


F 


22 X 10^^ 


34 


6.8 


G 


5.8 X 10^^ 


11 


2.2 


H 


19 X 10^^ 


4.0 


0.8 



nent nonthermal contribution. Indeed H04 show that the 
nonthermal fraction of the radio luminosity at 5 GHz is 
~ 0.7. If we adopt Model B with a burst SFR instead 
of Model A, the thermal component decreases, because 
the mean age of the stars is older. In this case also, how- 
ever, the nonthermal radiation of Model B falls well be- 
low the data presented in Figure [3 As long as we adopt 
the "standard" nonthermal component (i.e., Models A- 
D), we cannot explain the prominent nonthermal contri- 
bution. Therefore, we exclude Models A-D, and adopt the 
parameterized model for the nonthermal component (i.e., 
Models E-H). 

In the parameterized model of the nonthermal com- 
ponent, IntT-ai is treated as a parameter to be determined 
from observations fEa. I10|). Hence, we estimate ^ntTnt by 
using the observational data of SBS 0335— 052 taken by 
H04 and Dale et al. I|2001|l . We adopt the same initial 
conditions and parameters as above, and assume an age 
oi t — ^ Myr; we then search for the value of intent in 
Models E-H which minimizes for the six data points. 
The best-fit ZntTnt are listed in TableQlwith and (re- 
duced x^; the number of freedom is five), and the best-fit 
SEDs for Models E-H are shown in Figures Et-f, respec- 
tively. 

The burst models require higher l-ntT-at- In the burst 
models, all the stars form at the beginning, so that the 
mean stellar age is older. Thus, the thermal emission of 
the burst models is fainter than that of the continuous star 
formation models. This is why we require a higher non- 
thermal time-luminosity integral of a SN II in the burst 
SFR models (Models F and H) than in the continuous SFR 
models (E and G) to explain the observed radio emission. 
With the nonthermal component comparable to or larger 
than the thermal component, the points at 4.86, 8.46, and 
22.5 GHz are well reproduced. The spectral slope of the 
nonthermal component is also consistent with the obser- 
vations. Thus nonthermal radio emission comparable to 
or stronger than the thermal radio emission is required 
even in this young (< 5 Myr) galaxy. In general, the fits 
are satisfactory for Models E-H except for the point at 
14.9 GHz. The fractions of nonthermal emission are 0.44, 
0.79, 0.37, and 0.79 at 5 GHz in Models E, F, G, and H, 
respectively. The fractions in Models F and H are more 
consistent with the value of 0.7 at 5 GHz inferred from 
the spectral decomposition by H04. 

The nonthermal feature of the radio spectrum is in- 
consistent with an age oi t < 3.5 Myr, because of the 



time required for the onset of SNe. On the other hand, as 
long as we adopt an age between 4 Myr and 7 Myr, the 
difference in the derived Znt^nt is within a factor of ~ 2 
(i.e., within the uncertainty among Models E-H.) Thus, 
we conclude that Znt^nt derived here is robust against the 
age uncertainty. 

The nonthermal radiative energy emitted aA, v = 5 GHz 
over the entire lifetime of a SNR is predicted to be 5.8- 
22 X 10^2 W Hz-i yr (TableEJ. This is more than an order 
of magnitude larger than the value used in the "standard" 
model, in which the nonthermal radio energy emitted by 
a SNR is estimated based on the observed E — Z) (surface 
luminosity vs. radius) relation and the theoretical adia- 
batic lifetime. That assumption leads to the total radio 
energy per unit frequency as (assuming a = —0.5) 



1 cm 



El,i,.) = 1.1 X lO^^S^Y^" 

^(SGH^) ' 
while our new result suggests it to be 



-2/17 



i?P (i.) = (1.8-7.1) X 10 



30 



5 GHz 



-0.5 



J Hz" 



(28) 



(29) 



Integrating those quantities from 0.1 to 100 GHz, we ob- 
tain the total radio energy per SNR as 



(30) 



(31) 



£"(0.1-100 GHz) = 4.7 x lO^^E^^Y^" 

V 1 cm / 
for the "standard" case, and 

£P(0.1 - 100 GHz) = (8-31) X 10^^° J 

for the "parameterized" case. 

3.2. IZwlS 



We now apply our model to another "template" of nearby 
young BCDs, IZwl8, whose nebular oxygen abundance 
is low (12-Klog(0/H)=7.2, Skillman & Kennicutt 19®. 
The star-forming region of I Zw 18 is more diffuse than 
that of SBS 0335-052 (e.g.. Hunt et al. l:j()03all . Thus, the 
different gas density of IZwl8 from SBS 0335— 052 may 
provide us with independent information on the density 
dependence of radio emission. 

For the initial density and the radius, we assume riuo = 
100 cm~'^ and rsF = 100 pc, consistently with optical 
spectra and images (Skillman & Kennicutt Hunt 
et al. I2003at Hirashita & Hunt 200^. The resulting gas 
mass is Mgas = 1.4 x lO"^ Mq with tg = 4.4 x lO*' yr; the 
gas mass is comparable to the observationally determined 
values by van Zee et al. l(TU^ (2.6 x lO'^ Mq for their Hi- 
A component) and Lequeux & Viallefond H1980|l (2.1 x 
10^ Mq for their component 1). 

We adopt an age of 10-15 Myr for the dominant burst 
in I Zw 18. This is consistent with the results of our model 
of infrared dust emission (Hirashita & Hunt l^niM|l . which 
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Fig. 2. Radio SEDs at i = 5 Myr calculated by (a) Model A, (b) Model B, (c) Model E, (d) Model F, (e) Model G, 
and (f) Model H for the same initial condition as Figure ^ For Model E, we adopt the best fit parameter of ^nt''nt- 
The dotted and dashed lines represent the thermal and nonthermal components, respectively, and the solid line shows 
the sum of those two components. The points are observational data taken from Hunt et al. H2(J(J4|I (squares) and Dale 
et al. <j2001fl (cross). The reduced value (xl) is also shown in each panel. 
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gives an age of 10-15 Myr and consistent with the ob- 
served metalhcity and dust abundance. Takeuchi et al. 
iPMjl show that the dust mass derived by Cannon et al. 
H2()()2|l is consistent with an age of 10-30 Myr, which is also 
consistent with the observational upper limits of infrared 
flux. Recchi et al. H2004|l (and references therein) argue 
that stars older than 0.5-1 Gyr, if any, do not produce a 
significant contribution to the metal budget of I Zw 18, but 
intermediate age stars with an age of a few hundred Myr 
may be required from studies of the chemical abundances. 
An analysis of the color-magnitude diagram by Aloisi et 
al. H1999fl also suggests the existence of intermediate-age 
populations. Hunt et al. I|2003b|l derive an age of < 500 
Myr for the oldest stellar population (but see Ostlin 20QQ) 
and ^15 Myr for the youngest one in the main body. 

First, we evaluate a continuous SFR. The time evolu- 
tion of n, riH, EM, and SFR is shown in Figures |nt~d. 
In the first 1.4 Myr, the expansion speed of the ionizing 
front is faster than the sound speed because of small gas 
density. Thus, the dynamical response of the system be- 
gins later compared to SBS 0335— 052. Moreover, because 
of the low density, the emission measure remains smaller 
than in SBS 0335—052. The star formation stops at t ~ 7.9 
Myr because the entire star-forming region is ionized. The 
total stellar mass formed is 2.4 x 10^ Mq. 

In the burst SFR, the SFR is constant up to t^yn — 4.4 
Myr, a little shorter than the above duration (7.9 Myr). 
Thus, the expansion of the ionized region stops earlier 
and nn remains higher. However, the decline of rj occurs 
earlier because of the shorter duration, and the resulting 
EM is roughly similar to that for the continuous SFR. As a 
result, there is little difference in the evolution of emission 
measure between the continuous SFR and the burst SFR 
in a low-density regime. 

We also calculate radio SEDs with various models 
listed in Table As mentioned above, we adopt an age 
of 10-20 Myr. Because the (continuous mode of) star for- 
mation stops at t = 7.9 Myr, the radio emission decreases 
between t = 10 and 20 Myr. At 10 Myr, the thermal com- 
ponent is too high to be consistent with the observations 
reported by Hunt et al. H2()()5|l . The age should be older 
than 12 Myr to be consistent with the data. At 15 Myr, the 
thermal component is well below the observations. Thus, 
we examine the two ages, t — 12 Myr and 15 Myr, as high 
and low limiting luminosities for the thermal component. 
The emission measures at t = 12 Myr and 15 Myr are 
2.4 X 10^ pc cm~^ and 1.4 x 10^ pc cm~^, respectively. 
Those are larger than the EM estimated by Hunt et al. 
H2005|l (EM ~ lO'' pc cm"*^), but smaller than their ex- 
treme case (EM ~ 10^ pc cm~^). If we adopt an age of 
^ 20 Myr, the emission measure becomes too small to be 
consistent with observations. 

Now we examine the radio SED at t — 12 Myr, but 
consider only the p- models for the nonthermal component, 
since the s-models are rejected after the investigation of 
SBS 0335— 052. In fact, the s-models are also unaccept- 
able for IZwl8. As in Sect. 13.11 the best-fit nonthermal 
time-integrated luminosity (fluence) per SN II at 5 GHz 



Table 3. Best Fit Parameters for the Nonthermal 
Component (I Zw 18 at < = 12 Myr) 



Model /ntrnt(5 GHz) 
(W Hz-i yr) 



X 



E 
F 
G 
H 



1.1 X 10^ 
2.5 X 10^^ 
1.1 X 10^^ 
2.5 X 10^^ 



8.2 1.6 

7.5 1.5 
8.2 1.6 

7.6 1.5 



Table 4. Best Fit Parameters for the Nonthermal 
Component (IZwl8 at < = 15 Myr) 



Model /„trnt(5 GHz) ^ 
(W Hz"' yr) 



E 
F 
G 
H 



2.4 X 10^ 
3.0 X 10^ 
2.4 X 10^ 
3.0 X 10^ 



7.7 
8.5 
7.7 
8.6 



1.5 
1.7 
1.5 
1.7 



(IntTnt) is sought by minimizing for the six data points 
adopted from Hunt et al. f2005') and references therein. 
The data are shown in Figs. 01 and O the lower three data 
points are from the VLA (Hunt et al. 20[S|), while the 
upper three are from single-dish radio telescopes (Klein 
et al. 1991) . Although the discrepancy between those two 
data sets may arise from the sensitivity for the diffuse 
component, the discrepancy is only ~ 40% at ~ 5 GHz. 
Thus, the following conclusions are not changed by the 
uncertainty in the diffuse component. 

We include all the six data points to obtain the best-fit 
values for IntTnt, which are shown in Table |21 for t — 12 
Myr, and in Table 01 for t = 15 Myr. Applying the values 
reported in those tables, we examine the radio SEDs. The 
results are shown in Figs. 01 and for t = 12 and 15 Myr, 
respectively. Those figures confirm that the nonthermal 
luminosity derived from the parameters in Tables O and 
01 is consistent with the data points. Moreover, the spec- 
tral index at i/ ^ 5 GHz is reproduced by introducing a 
significant contribution from the nonthermal component. 

The spectral decomposition by Hunt et al. (2005 ) gives 
a nonthermal fraction of 0.59 at 5 GHz, corresponding to 
a thermal fiux of 0.54 mJy. This is approximately realized 
for Models F and H shown in Panels b) and d) of Fig. 
01 and for Models E and G shown in Panels a) and c) of 
Fig. |S1 However, other solutions with a thermal compo- 
nent of 0.3-1 mJy are also permitted. The younger age of 
12 Myr gives 0.20 Mq yr~^ for the averaged SFR, while 
15 Myr gives 0.16 Mq yr~^ . The massive SFR is estimated 
as SFR>5 Mq ^ 0.036 and 0.029 Mq yr-\ respectively. 
Those values are 1.5-2 times larger than the estimate of 
Hunt et al. H2005|l . If we assume a smaller SFR according 
to Hunt et al. H2005|l . ^nt^nt required to fit the observa- 
tional SED becomes 1.5-2 times smaller. 

Cannon et al. H2005|l also provide fluxes of I Zw 18 at 
V — lA GHz and v = 8.4 GHz; their fluxes are lower than 
the values given by Hunt et al. H2005|l . but it is possi- 
ble that the discrepancy comes from the sensitivity to the 
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diffuse component. However, they also derive a spectral in- 
dex that indicates a prominent nonthermal contribution. 
Consequently, our conclusion that a strong contribution 
from the nonthermal component is required should be ro- 
bust (and ^ntTnt derived by us should be relatively secure). 

We have derived values of ZntTnt for IZwl8 between 
1.1 X 10^2 W Hz-i yr and 2.5 x lO^^ W Hz"i yr at 5 GHz. 
Since the radio spectrum is dominated by the nonthermal 
component already at 15 Myr, those derived values should 
be robust for older ages, because the thermal component 
steeply declines thereafter. On the other hand, if we adopt 
an age smaller than 10 Myr, the fraction of the thermal 
component is too large to be consistent with observa- 
tions. Thus, even though there is uncertainty in the age of 
I Zw 18, we can safely conclude that /ntTnt of I Zw 18 is two 
to five times smaller than that of SBS 0335—052. The typ- 
ical densities of ionized regions are ~ 10"^ cm~'^ and ^ 10^ 
cm-3 in SBS 0335-052 and I Zw 18, respectively. We pro- 
pose in Sect. 21a scaling relation ZntTnt oc 7^0-35~o.73^ Since 
^jjjTnt = 5.8-22 X 10^2 W Hz~^ yr has been obtained for 
SBS 0335—052, the above density scaling relation predicts 
that IntTnt ^ 1.1-9.9 X 10^2 W Hz"! yr. The above values 
derived for I Zw 18 arc in this range. 



For IZwl8, we obtain the following value for the non- 
thermal radiative energy over the entire lifetime of a SNR: 

= 3.4-9.4 X lO^s (^^) J Hz-' . (32) 

As for SBS 0335— 052, this is larger than the "standard" 
estimate in Eq. 128|l . The total radio energy per SNR in 
I Zw 18 becomes 

£P(0.1-100 GHz) = 1.5-6.0 X 10^^° J . (33) 

4. Synchrotron emission and magnetic fields 

Through comparison with observations, we have obtained 
the radio energy radiated by the nonthermal synchrotron 
component. In this section, we relate the observationally- 
derived values to the physical quantities governing the 
synchrotron radiation. 

4.1. General estimate 

The energy spectrum of energetic electrons radiating syn- 
chrotron radiation is assumed to follow a power law as 

N{E) dE = CE-P dE , (34) 
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Fig. 4. Radio SED at t = 12 Myr calculated for IZwl8 by using Models E, F, G, and H (panels a), b), c), and d), 
respectively). The dotted and dashed lines represent the thermal and nonthermal components, respectively, and the 
solid line shows the sum of those two components. The squares are observational data taken from Hunt et al. H2[)()5|l 
and references therein. The reduced value (xl) is also shown in each panel. 



where N{E) dE is the number density of the energetic 
electrons with energies between E and E + dE, and C 
is the normalizing constant. Under this energy spectrum, 
the energy of the energetic electrons per unit volume, Ec, 
is calculated as 

Ec = / EN{E) dE , (35) 

where Emax and i?min are the maximum and minimum 
energies of the energetic electrons. Combining equations 
H34|) and (|35|l , we can express the normalizing constant as 



£c(2-p) 



C 



E, 



2-p 



E: 



2-p 



(Py^2) 
(p-2) 



(36) 



ln(i?max/£^min) 

The emissivity of the synchrotron radiation, P{v), 
is written as (Ginzburg & Syrovatskii 119651 Rybicki & 
LightmanCnZni) 

P{i^) = 1.4 X 10'2^a(p)CB(''+i)/2 



6.3 X 1018 Hz 



-(p-l)/2 



(37) 



where a(p) is a coefficient which depends on the exponent 
p in the energy spectrum defined in equation 134|l . C is 
the normalizing constant in equation H34() . and B is the 
magnetic field strength in units of gauss. 

The spectral index observed in individual SNRs is 
roughly a ~ -(p - l)/2 ~ -0.5 (Clark & Caswell ITWHIl : 
i.e., p ^ 2. Since the frequency of the power-law syn- 
chrotron radiation ranges at least from MHz to GHz, we 
can assume that ln(ii^max/-E'i„in) ~ 10. In any case, be- 
cause of the logarithmic dependence, the results are not 
sensitive to the assumed £^max/-6'min- 

The total nonthermal radio luminosity of a region can 
be calculated by integrating P(z/) over all the volume V 
as 



r theory 
^nt 



/ 47rP(l/)dF, 

Jv 



(38) 
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where LJJ""'^'^ is the theoreticaUy estimated nonthermal 
luminosity of the entire region. Assuming homogeneity of 
the region and p = 2, and with a(2) ~ 0.103 (Ginzburg & 
Svrovatskii 11965(1 we obtain 



retical expression for IntTnt ai v = 5 GHz as 



^ theory /X _ 



(zy) = 2.0 X 10 



19 



n 

V 1^ 



1052 erg J V 10 AiG 



B 



1.5 



W Hz-1 , 



,5 GHz/ 

where is the total energy in the volume: 

= [ e,dV. 
Jv 



(39) 



(40) 



We assume that a fraction /e of the input SN energy is 
converted to the energy of accelerated energetic electrons: 



£c — fcEsTSlJTc , 



(41) 



where Te is the lifetime of energetic electrons. Equating 
L°t p (Eq. and L*t™'^ (Eq. EHI), we obtain the theo- 



ZntTnt = 2.0 X 10 



19 



ln(i?max/^^min)\ ^ / fcE; 



105 yr y V 10 mG 



B 



10 J \ 10-17 erg 

1.5 



(42) 



It would be reasonable to assume that r^t — Tq. However, 
we retain both quantities since the observational con- 
straint has been applied to IntTnt- 

4.2. Comparison with SBS 0335- 052 

In Sect. Ofor SBS 0335-052 we derived the nonther- 
mal radio energy emitted at 5 GHz per SN as ?nt''nt — 
(6-22) X 10^2 W Hz~^ yr. Comparing this value with 
the theoretical expression given in Eq. (|42|l . we obtain 
(/e/10-4)(Te/105 yr)(B/10 AiG)i-5 ^ (3-11) X 10^ where 
we have assumed i?sN — 10^^ erg and ln(i?max/-£'min) = 
10. 

Either B > 10 /iG or /e > lO^^^ or both. We first in- 
vestigate the possibility that the magnetic fields could be 
amplified by turbulence generated from SNe as proposed 
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Fig. 6. Time evolution of an "active" radio SED. The same parameters as used to model SBS 0335—052 are adopted 
as a representative case for the "active" mode. Panels a), b), c), and d) represent the SEDs at 3.5, 5, 10, and 15, 
respectively. The dotted and dashed lines represent the thermal and nonthermal components, respectively, and the 
solid line shows the sum of those two components. 



by Balsara et al. f 20041. We define the ratio of thermal en- 
ergy to magnetic energy, (3 = [ink^T / 2) / {B'^ / Sir) , where 
n is the particle number density, k-Q is the Boltzmann con- 
stant, and T is the gas temperature. In SBS 0335—052, we 
can assume n ^ 10^ cm"^ and T 2 x 10^ K (Sect. ETJ. 
Then, the magnetic field strength becomes 320/3^^/^ ^G. 
The magnetic field could be amplified in the star-forming 
region until the magnetic pressure becomes comparable 
to the thermal pressure (otherwise, the magnetic fields 
would not be confined within the region). Thus, (3 > 1 
would be a reasonable condition. If we conservatively as- 
sume that (3 ~ 10, we obtain the magnetic field strength 
in SBS 0335-052 as 100 /.tG. The production efficiency of 
energetic electrons is not well known, but it could be as- 
sumed to be of an order of 1% or more (e.g., Zirakashvili 
& V61k'2006). Here we conservatively assume /c ^ 0.01. 
With the quantities estimated in this paragraph, we pre- 
dict Tc - (6 - 35) X lO'' yr. 

This timescale can be interpreted in two ways. One is 
the radiative lifetime of a SNR. The adiabatic lifetime of 



a SNR is estimated as Tad ^ 10'^ yr if we assume nn ~ 10^ 
cm~'^ (Eq. inj. This is much shorter than the above To. 
However, the SNR continues to expand after the cooling 
due to the momentum conservation. If we adopt the ex- 
pansion law of ShuU H1980fl . it finally reaches the sound 
speed of the ionized medium and the shock disappears af- 
ter ~ 10^ yr. Therefore, the production of energetic elec- 
trons in shocked gas may continue on a timescale of ~ 10^ 
yr, which is roughly consistent with Te estimated above. 
From the expansion law, we can derive Te oc 

The other interpretation of Te is that it reflects the life- 
time of an energetic electron. The energy loss of energetic 
electrons should be considered at frequencies larger than 
v-m estimated as 



3.1 X 10^3 



(43) 



where B±^ is the component of the magnetic field per- 
pendicular to the velocity of the electron, and the time 
t is measured from the production of the energetic elec- 
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trons (Ginzburg & Svrovatskii I1965f) . We define Tg as the 
timescale on which the energy loss of electrons contribut- 
ing to the radiation at a frequency v becomes significant. 
Then, Tg can be estimated by inserting z/„i = v and t = Ts 
in Eq. as 



8.0 X lO*" 



10 



-3/2 



GHz 



-1/2 



(44) 



If we adopt B± - 100 ^G for SBS 0335-052, wc obtain 
Ts ^ 2 X 10^ yr at ^ 10 GHz. This is consistent with the 
above Tnt . In the following, we adopt an approximation of 



B. 



4.3. Scaling ofl^tTnt 

We investigate both cases for Tnt. The first is based 
on the radiative lifetime of a SNR, i.e., Tnt = ~ 
10^(71/10"^ cm^'^)^^/^ yr, and the second on the lifetime of 
energetic electrons, i.e., Tnt = Ts given in Eq. (|^ . By us- 
ing Eq. H42(l , these scaling relations result in the following 
expression for IntTnt- 



^nt^Vit 



21 f fcEs 



2.0 X 10^1 ( 

V 1049 erg 

B \^-^ 
X I I WHz-i 



) (l03 cm-s) 
yr 



-2/5 



10 AiG, 

,23 ( f'^EsN 



(45) 



3.6 X 10^ 



V 10*9 erg 

if TSN = Ts, 



WHz- 



yr 



where we have assumed that ln(£'niax/-E'„iin) — 10. 

It is generally difficult to observe the magnetic field 
strength in BCDs. If we use /3 (Sect. 14.21) . we can ex- 
press the magnetic field strength in terms of the gas den- 
sity as B ~ 320 /iG (n/lO^ cm-3)^/2^-i/2^ where we 
have used T = 2xlO^K.A similar observational scal- 
ing {B oc n°'^®) was observationally found by Niklas & 
Beck (^97). Inserting this expression into the first case 
in Eq. H45() . we obtain a scaling relation as ^nt^Vit = 8.1 x 
1023(/^£;sn/1G49 erg)(n/103 cm^^^o.ss^-o.rs ^ jjz"! yr. 
If we assume that /ci?SN does not change among galax- 
ies and that ^ntTVit scales only with density, we obtain 
tnt Tnt ^ (^ntTnt)sBs("/10^ Cm 3)0-35^ where (ZntTnt)sBS is 
the value obtained for SBS 0335—052 as reported in Table 

m 

In the second case in Eq. I|45(l , there is no dependence 
on magnetic fields and gas density, if /c£^sn does not de- 
pend on those quantities. Thus, we also consider the case 
where IntTnt is independent of n. 

Finally we propose yet another scaling of the magnetic 
field strength, assuming that the magnetic field is am- 
plified by SNe, as in the turbulence-driven amplification 
mechanism proposed by Balsara et al. H2004fl . Since the 
SN rate is proportional to the SFR, the energy density 
of the magnetic fields could scale with the SFR density 
('- SFR/r^). Equation imphes that the SFR density 
scales as n^^^ (this scaling also expresses a Schmidt law; 



Schmidt T^). Thus, we obtain B oc n^/^. With this scal- 
ing relation, and normalizing it to SBS 0335— 052, we ob- 
tain ^ntTnt — (^ntTnt)sBs("/10'^ ciR^^)^ '^^ for thc first case 
in Eq. (|45(l . Since the energy density of energetic electrons 
is proportional to the SN rate, this scaling indicates that 
the magnetic energy density scales with the energy density 
of energetic electrons. 

In summary, the possible scaling of ^ntTnt is summa- 
rized by the following three cases: 

(^ntTnt)sBS Scahug a, 

IntTnt = < (?ntTnt)sBs(n/103 cm-3)0-35 Scaling b, (46) 
(^ntTnt)sBs("/10' cm-3)0.73 Scaling c. 

Indeed, it has recently been found that the S-_D (radio 
surface brightness to diameter) relation depends on the 
ambient density of interstellar medium (Arbutina et al. 

Arbutina & UrosevicESS) : SNRs with higher ambi- 
ent density tend to emit more radio energy. Thus, the pos- 
itive correlation between density and IntTnt is supported. 
We should note that thermal emission may contribute to 
the radio emission from SNe especially in dense environ- 
ments (Urosevic & Pannuti l2005|l . Thus, our derived val- 
ues for ZntTnt may include some contribution from thermal 
emission. The decomposition of thermal and nonthermal 
contributions is left for future work. 



4.4. Comparison with IZwlS 

The typical number density of IZwlS is n c± 100 cm^'^ 
(Sect.E3l- Comparing Table|a(model for SBS 0335-052) 
with Tables El and m (model for IZwlS at 12 and 15Myr, 
respectively), we find that the best model estimate for 
(^ntTnt)iZw/(^ntTnt)sBS — 0.13 - 0.33, whcrc the subscript 
"IZw" indicates the value for IZwlS. Scaling a would 
imply no density dependence, so that we would expect 
(^ntTnt)iZw — (^ntTnt)sBs; this is highly inconsistent with 
our models. If, instead, we use Scaling b (oc nP'^^), we 
would derive (/ntTnt)iZw/('ntTnt)sBS — 0.45, stiU a bit 
larger than our model predictions. Scaling c (oc nP''^^) 
would give (ZntT„t)izw/(^ntTnt)sBS - 0.19, which is the 
most consistent with our results, although Scaling b can- 
not be rejected with certainty. 

Interestingly, the nonthermal luminosity observed in 
discrete SNRs in ten external galaxies appears to be pos- 
itively correlated with density (Hunt & Revnolds 12006(1 . 
in a way that is consistent with our results. If the time- 
luminosity integral varies as n"'^, as in the seemingly 
most likely dependence, then the observed correlation of 
Int OC n^'^ implies that the radiative lifetime of a remnant 
should vary as The slope of the Int — n relation 

is rather uncertain (Hunt & Revnolds 12006(1 . but cannot 
be less than unity. With Int esc n, the remnant lifetime 
would vary as n~^'^. These values encompass the expan- 
sion lifetime dependence Tc oc n~^/^ and the adiabatic one 



To oc n 
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5. The time evolution of the radio continuum 



calculated by the following equation: 



Cannon & Skillman ( 2004) propose the radio spectral in- 
dex as an age indicator. The radio spectral index is de- 
fined by fitting the radio spectrum with a functional form 
of v"^ . In their paper, the observational spectral index is 
determined with the data at v = 1.4, 4.9, and 8.5 GHz. 
They argue that af > with ages ^ 1 Myr because of the 
free-free absorption. Then the optically-thin thermal spec- 
tral index (af = —0.1) appears within a typical lifetime of 
H II regions (~ 10 Myr) and finally the spectrum steepens 
because of the increasing contribution from a synchrotron 
component. 

However, we have shown in Sects. IXTl and 15^ that the 
gas density is also important for the shape of the radio 
spectrum. The dense and compact star-forming region in 
SBS 0335—052 shows a strong burst of star-formation with 
strong free-free absorption. We have dubbed this char- 
acteristic "active". On the other hand, the diffuse star- 
forming region in I Zw 18 has more mild star-formation 
with much less absorption. We have called this mode "pas- 
sive" . Since the radio free-free absorption affects the radio 
spectral index, we can distinguish these two classes of star 
formation, "active" and "passive" , through the radio spec- 
tral index. This also means that the age is not the only 
factor that determines the spectral index; the gas density 
is also important in its determination. 

Here we examine the time evolution of radio spectra, 
adopting the same parameters as used for SBS 0335— 052 
and IZwlS, considering them as representative cases for 
"active" and "passive", respectively. Model G is adopted 
here, but the qualitative behavior of the radio spectral 
index is not much affected by the specific details of the 
adopted models. We use IntTnt estimated for Model G in 
Tables El and El for the active and passive cases, respec- 
tively. In Figs. El and |3 we present the radio SEDs cal- 
culated by the models for SBS 0335-052 and IZwl8, re- 
spectively. The snapshots at t — 3.5 Myr (soon after the 
first explosion of SNe at ~ 3 Myr), 5 Myr, 10 Myr, and 15 
Myr are shown. We observe that the free-free absorption 
around v = 1 GHz is prominent in the active case even at 

5 Myr. On the contrary, there is little absorption even 
at t ~ 3.5 Myr in the passive case. This indicates that 
the absorption feature is not a simple indicator of age as 
suggested by Cannon & Skillman f20Q4') . The density of 
the star-forming region is also important for the spectral 
index around ly ^ 1 GHz. 

By comparing Figs. El and we also see that the 
nonthermal component becomes comparable to the ther- 
mal luminosity earlier in the active case than in the pas- 
sive case. This mainly comes from the difference in time- 
integrated radio luminosity of a SN (Znt'^nt)- Thus, we sug- 
gest that a nonthermal-dominated radio SED is observed 
more often in active BCDs than in passive BCDs. This 
is true if IntTnt correlates with gas density. Observations 
compiled by Hunt & Reynolds l)2006|l support this view. 

Finally we show the time evolution of radio spectral 
index af. The spectral index defined bX v — vi and 1^2 is 



m{vi, V2) 



\og[L{v^)/L{v2)] 

log(j/i/l/2) 



(47) 



We put vi — 1.5 GHz and V2 = ^ GHz to be consistent 
with the frequency range often used to determined the 
spectral index.Figure El shows the results. The solid and 
dotted lines represent the results of the "active" and "pas- 
sive" models whose spectral shapes are shown in Figs. El 
and O respectively. The spectral index depends on age, 
and it approaches the nonthermal value (—0.4) as the age 
increases. This confirms the trend suggested by Cannon & 
Skillman H2U04|I . More importantly though, the spectral 
index is not a simple function of age but rather depends 
strongly on the density of the star- forming region (i.e., 
"active" and "passive" modes). Indeed the spectral index 
is positive even at t = a few Myr in the active mode, while 
it is always negative in the passive mode. 

Another important consideration is that the radio 
spectrum can be flat (i.e., at ~ 0) even when the con- 
tribution from the nonthermal component is significant. 
Indeed, at t = 5 Myr, the nonthermal and thermal com- 
ponents are comparable in the active mode (Fig. EJ, and 
the spectral slope at high frequency indicates the presence 
of the nonthermal contribution. However, if we define the 
spectral index by using v — 1.5 GHz and 5 GHz, the 
spectral index becomes ~ 0. This spectral index could be 
misinterpreted as thermal. Thus, in order to avoid such 
a misinterpretation, it is important to derive a spectral 
index at ^ 5 GHz. 



5.1. Implications for primeval galaxies 

High-z primeval galaxies tend to have a high density; 
they are formed and evolve in deep gravitational poten- 
tials (e.g., Padmanabhan [n)93fl . Therefore, the gas in such 
primeval objects may have high pressure, and the star- 
forming regions may be expected to mimic the active mode 
(Hirashita & Hunt l^nUljl . Indeed, the number density of 
gas in high-z galaxies is estimated as ^ 10'^ cm""^ (Norman 
& Spaans lTWTjl . and is similar to that of SBS 0335-052, 
typical of the active class. Therefore, the nonthermal com- 
ponent could dominate on short timescales after the onset 
of star formation while the radio emission is optically thick 
at <^ 1 GHz. The dense environment could also aid am- 
plification of magnetic fields as proposed by the scaling 
relations in Sect. 14.31 

This paper also has implications for the production 
of high-energy electrons. It would be reasonable to con- 
sider that the production occurs on an expansion timescale 
SNRs, since the shocks associated with SNRs can be re- 
sponsible for the acceleration of electrons (e.g., Axford 
11981(1 . The lifetime of nonthermal radiation per SN is de- 
termined by the minimum of two timescales: the lifetime 
of the shock and the synchrotron loss timescale. 
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Fig. 7. Same as Fig. but for the "passive" mode, with the IZwl8 model as a representative case. 



6. Conclusions 



In order to investigate the time evolution of the radio 
SEDs of metal-poor star- forming galaxies, we have con- 
structed a radio SED model by treating the thermal and 
nonthermal components consistently with the star forma- 
tion history. In particular, the duration and luminosity 
of the nonthermal radio emission from supernova rem- 
nants (SNRs) has been constrained by using the obser- 
vational radio SEDs of SBS 0335-052 and IZwl8, both 
of which are reasonable proxies for young galaxies in the 
nearby universe. In SBS 0335— 052, the typical radio en- 
ergy emitted per SNR over its lifetime is estimated to be 
£{0.1 - 100 GHz) - 8-32 x lO"*^ erg. In IZwl8, another 
representative of metal-poor star-forming galaxies, we find 
that £:(0.1 - 100 GHz) ~ 2-6 x lO'^'^ erg. Both estimates 
are significantly larger than previous estimates based on 
the E — _D relation. These values of the two "template" 
galaxies can be simultaneously explained by a simple den- 
sity scaling relation of i? oc 77,0-35-0.73^ indicating that the 
magnetic pressure scales with the gas pressure, or that the 
magnetic fields are amplified by SNe, or both. 



We have also predicted the time dependence of the ra- 
dio spectral index and of the radio SED itself, for both the 
active (SBS 0335— 052) and passive (IZwlS) cases. These 
models enable us to roughly age date and classify radio 
spectra of star- forming galaxies into active/passive classes. 
Since the radio emission around i/ ^ 1 GHz can be affected 
by the free-free absorption especially in the active class, 
the spectral index defined at ^ 5 GHz should be used 
to estimate the contribution from the nonthermal compo- 
nent. On the other hand, the free-free absorption feature 
around ~ 1 GHz could be used to select the active class. 
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Fig. 8. Time evolution of the radio spectral index defined 
at 1.5 GHz and 5 GHz [a(1.5 GHz, 5 GHz)]. The sohd and 
dotted lines show the results of the "active" and "passive" 
models (same as Figs. El and [7| respectively). 
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Astrophysics Data System Abstract Service (ADS). 

Appendix A: Analysis of expanding H ii regions 

We derive some approximate analytic solutions for ex- 
panding ionized regions based on Eqs. (|23(l and H24|l . By 
combining these two equations, we obtain 



3/4 



(A.l) 



Since we estimate rjo by the Stromgren radius under the 
current stellar luminosity and the gas density of the neu- 
tral medium outside the ionized region, we obtain 



f io = A 
where 

A = 



2/3 



47r(1.08)a(2) 



1/3 



(A.2) 



(A.3) 



The factor 1.08 comes from the correction for ionized he- 
hum (Eq. O). By combining Eqs. HA.l|) - (jA.3|l . we obtain 



drj 
dF 



Solving this equation, we obtain 



(A.4) 



(A.5) 



where we have assumed that ri(t 0) = 0. If we assume 
that the SFR is constant and that the death of stars is neg- 
ligible, iVion is almost proportional to t, and we can pre- 
dict that n oc i^/^. This relation, combined with Eqs. (|23|l 
{pii oc (dri/dt)2), indicates that pu oc t'^/"^ . The emission 



measure evolves as EM oc pjjrj oc t~^^'' . In spite of this 
rough derivation, those relations explain the behaviors de- 
picted in Fig. n In the burst SFR (Fig. 0, the emission 
measure is particularly deviates from the predicted power- 
law evolution. This is because the first increase of n is due 
to increase of the Stromgren radius itself, not due to dy- 
namical expansion. 
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